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We provide a new approach to open quantum systems which is based on the Feshbach projection 
method in an appropriate Hilbert space. Instead of looking for a master equation for the dynamical 
map acting in the space of density operators we provide the corresponding equation for the evolution 
in the Hilbert space of the system. Its solution enables one to construct a legitimate quantum 
evolution (completely positive and trace preserving). The big advantage of this approach is the 
ability of performing a consistent approximation which creates notorious problems in the standard 
Nakajima-Zwanzig approach. The new scheme is illustrated by well known spin-boson model beyond 
rotating wave approximation. 

PACS numbers: 03.65.Yz, 03.65.Ta, 42.50.Lc 



a* 



(N 
> 
oo 

(N 

\o 

O 

m 



x 



Introduction. - The description of a quantum system 
interacting with its environment is of fundamental impor- 
tance for quantum physics and defines the central objec- 
tive of the theory of open quantum systems (HQ- During 
the last few years there has been an increasing interest in 
open quantum systems in connection to the growing in- 
terest in controlling quantum systems and applications in 
modern quantum technologies such as quantum commu- 
nication, cryptography, computation and ever growing 
number of applications. In practice, this theory is usu- 
ally applied in the so-called Markovian or memoryless 
approximation. However, when strong coupling or long 
environmental relaxation times make memory effects im- 
portant for a realistic description of the dynamics one 
needs more refine approach and hence the general struc- 
ture of non-Markovian quantum evolution is a crucial 
issue. For the recent papers devoted to both theoreti- 
cal and experimental aspects of quantum evolution with 
memory see e.g. a collection of papers in Q and refer- 
ences therein. 

The standard approach to the dynamics of open sys- 
tem uses the Nakajima-Zwanzig projection operator tech- 
nique [U which shows that under fairly general condi- 
tions, the master equation for the reduced density matrix 
p(t) takes the form of the following non-local equation 



dt 



Pi 



K-t-uPu du 



(1) 



in which quantum memory effects arc taken into account 
through the introduction of the memory kernel /Q: this 
simply means that the rate of change of the state p t at 
time t depends on its history. An alternative and tech- 
nically much simpler scheme is the time-convolutionless 
projection operator technique [H, [H, 0] i n which one ob- 
tains a first-order differential equation for the reduced 
density matrix 



dt 



Pt = L t p t . 



(2) 



The advantage of the local approach consists in the fact 
that it yields an equation of motion for the relevant de- 
grees of freedom which is local in time and which is there- 
fore often much easier to deal with than the Nakajima- 
Zwanzig non-local master equation ([1]). 

It should be stressed that the structure of the memory 
kernel ICt is highly nontrivial and, therefore, the non- 
local master equation ([T]) is rather untractable. Note, 
that this equation is exact, i.e. in deriving ([T]) one does 
not use any specific approximation. Approximating ([1]) is 
a delicate issue. One often applies second order Born ap- 
proximation which considerably simplifies the structure 
of /C t . However, this approximation in general violates 
basic properties of the master equation like for example 
complete positivity or even positivity of pt @. Further 
simplification of (JT|) consists in various Markov approxi- 
mations which allow one to avoid memory effects. These 
approximations may also break the physics of the prob- 
lem. For example well known local Redfield equation 
@ again violates complete positivity 0, The prob- 
lem of a consistent Markov approximation was studied 
by Diimcke and Spohn One often tries to use phe- 
nomenological memory kernels. However, as was already 
observed by Barnett and Stenholm [l2j|, there is no sim- 
ple recipe how to construct ICt in order to preserve basic 
properties of quantum evolution [l3|, Q • 

The local approach based on ^ is much more popu- 
lar. Clearly, the mathematical structure of ([2]) is much 
more "friendly" than that of ([1}. Moreover, the local 
equation provides a straightforward generalization of the 
celebrated Markovian semigroup governed by the time- 
independent generator L possessing well known standard 
form [H [lj] 

Lp = -i[H, P] + \Y. V ^ + [ V ^PV1]) > (3) 

with the effective Hamiltonian H and arbitrary noise op- 
erators V a . Note, however, that local equation ((2]), con- 
trary to generalized master equation ([1]), is not derived 
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from any specific microscopic model. Therefore, it is a 
rather phenomenological equation which generalizes ([3]) 
for time-dependent situation. Again, there is no general 
recipe how to construct legitimate local generators. It 
should be also stressed that Markovian semigroup, being 
a special case of J2]) , is derived from (JT]) by applying quite 
sophisticated Markovian approximations like for example 
weak coupling or singular coupling limits [ToL ITl| . 

Recently, there are several proposals to define and de- 
scribe non- Markovian evolution (l7l - [20| . Rivas et. al. 
fl8| define non-Markovianity of the evolution by devia- 
tion from divisibility. A different approach is advocated 
by Breuer et al [ljj. It defines non-Markovian dynam- 
ics as a time evolution for the open system characterized 
by a temporary flow of information from the environment 
back into the system and manifests itself as an increase in 
the distinguishability of pairs of evolving quantum states. 
According to this definition the evolution is Markovian 

^ c§H'°* — "til 1 — f° r an initial states p and a and 
t > 0. Optimizing over p, a enables to construct suitable 
non-Markovianity measure (l9j (for comparing these two 
notions of non-Markovainity see e.g. [21H23| ). Other 
measures of non-Markovianity were recently proposed as 
well IH-il. 

A different approach to reduced dynamics of open 
systems is based on the solution of the Schrodingcr 
equation for the original composed system: one solves 
Ji\^t) = —iH\tyt) and then defines reduced dynamics 
Pi = tTE \i$t)(*&t\ by tracing out environmental degrees 
of freedom. It is clear that such solution gives us only a 
single trajectory in the space of density operators. Note, 
however, that if we start with \ty) = (g> \i/je) as an ini- 
tial state, then linearity of evolution enables us to con- 
struct a legitimate dynamical map A t (po) for any initial 
mixed state po of the system. In this case one easily finds 
that formally p t = A t (po) satisfies the local equation ([2]) 
with time dependent generator defined by L t := AtAj" 1 . 
By construction this generator is physically legitimate 
giving rise to complete positive and trace preserving dy- 
namical map At- Note, however, that this construction is 
slightly artificial since knowing A t the problem of reduced 
dynamics is completely solved even if the corresponding 
equation for At is not know. The prominent example of 
this approach is well known spin-boson model 

Gauge theory of quantum systems. - In this Letter we 
provide a new approach to the reduced dynamics of open 
quantum systems. Instead of applying the Nakajima- 
Zwanzig projection we apply the Fcshbach projection for- 
malism [29( to the Schrodingcr equation of the total sys- 
tem. This formalism was recently app lied in the con- 
text of open quantum system in [3CH32J . Our approach, 
however, is based not on the Schrodingcr equation for 
the vector state \^f t ) but on the Schrodinger-like equa- 
tion for the "amplitude" operator. Let as recall that if p 
is density operator in H then a Hilbert-Schmidt opera- 
tor k e C 2 (H) is called an amplitude of p if p = kkL 
Recall, that £ 2 (H) is equipped with a scalar product 
(k, rf) — tc^ri) and k e £ 2 (H) if the Hilbert-Schmidt 



norm ||k|| 2 = (k, k) is finite. Note that tip = 1 implies 
(k, k) = 1. If = a is an observable, then 

(k, a k) = tr (k^ok) = tr (ann^) = tr (ap) , (4) 

reproduces the standard formula for the expectation 
value of a in the state p. Amplitudes display a natu- 
ral gauge symmetry: a gauge transformation k — >• nU 
leaves p = kk^ invariant for any unitary operator U in 
T~i. The main idea of this paper is to analyze the dynam- 
ics of gauge invariant p in terms of its gauge dependent 
amplitudes k. This is well known trick in physics. Recall 
that for example in Maxwell theory it is much easier to 
analyze Maxwell equations not in terms of gauge invari- 
ant E and B fields but in terms of gauge dependent four 
potential A^. Suppose that p t satisfies von- Neumann 
equation 

idtpt = [H t ,pt] , (5) 

with time-dependent Hamiltonian H t . Its solution is 
given by p t = VtpV t , where the unitary operator V t 
solves the Schrodinger equation idtVt = H t Vt with the 
initial condition Vo = I. The corresponding amplitude n t 
satisfies 

id t K t = H t K t - K t G t , (6) 

where Gt is an arbitrary time-dependent Hermitian op- 
erator. One finds Kt = VtnW} , where the unitary oper- 
ator Wt solves the Schrodinger equation idtWt — GtWt- 
Hence n t n\ = Vtnn^ . If Gt = H t then n t satisfies the 
same von-Neumann equation as p t (such choice is used 
e.g. in (HI). Let us observe that if Kt solves ^ then a 
local (i.e. time-dependent) gauge transformation 

ti' t = K t Ut , (7) 

does not effect p t = n' t K 't ■ Note, that n' t solves the fol- 
lowing equation 

id tK ' t = Htn't - K ' t G' t , (8) 

where 

G' t = u}GtUt - lU^dtUt , (9) 

and hence Gt transforms like a non-abelian gauge field. 
It is, therefore, clear that taking U t such that idtUt = 
GtUt one kills G-term and arrives at the Schrodingcr-likc 
equation for an amplitude 

id t K t = H t K t . (10) 

As we shall see this simple choice leads to considerable 
simplification of the underlying structure of the gauge 
theory. Note that the Schrodinger-like equation (ff"Tj|) still 
allows for global (i.e. time independent) gauge transfor- 
mations K t —¥ K t U. 

Feshbach projection technique. - Consider a quantum 
system coupled to the environment living in %$ ®"He 
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and let H denote the total Hamiltonian of the composed 
system. Moreover, let p(g>f2 be an initial state of the 
composed system and denote by Kg ® ke an amplitude 
of the initial state, i.e. p = KgK s an< ^ ^ = k ek e . Wc 
introduce orthogonal projectors P a : £ 2 (Hg ®He) — > 
£ 2 (Hs®H E ) 

Pops® He = Hs® k e {k e ,Pe) , Px = l-P . (11) 

Let K t e C 2 (Hg ®'He) satisfy Schrodinger-like equation 

id t K t = HK t , Kq = Kg ® k e , (12) 

and (ksi Kg) = (ke, he) = 1- Standard Fcshbach projec- 
tion technique |29| leads to 

d t PoK t = -iHo P K t - iHoiP\Kt , (13) 
d t PiK t = -iH 1Q P Kt ~ iHuPiKt , (14) 

where we introduced a convenient notation = PiHPj . 
Solving for the irrelevant part P\K t 

PiK t = -i [ dse- l{t - s)Hll H 10 P K s , 
Jo 

one ends up with the following non-local equation for the 
relevant (system) part Po^t'- 

dtPo^t = -iHo P K t - / K t _ s P K s ds , (15) 
Jo 

with K t = H 01 e-' ltHll H w . Let Z t :U s -> H<? be defined 
by 

(Z t Kg) <S> k e = PoK t = PoU t (K S ® K E ) , (16) 

where f/t = e~ %m . Equation (fT5")) may be rewritten as 
the following equation for Z t 

dtZ t = -M eS Z t - f M t - s Z s ds, (17) 
Jo 

where the effective system Hamiltonian is defined by 
H ctf = tr E (H • I®fi) and M t = tr E (K 4 • I® 12). The 
above equation may be formally solved in terms of the 
Laplace transform Z p as follows Z p = (pl+iHes+Mp) -1 . 
Let us observe that (Z t Kg, Z t Kg) < (Kg, Kg) which shows 
that Z t Kg is no longer a legitimate amplitude for t > 0. 
It is clear since Z t describes the decay of Kg and hence 
|ZfKs|| is not conserved - it leaks out to the irrelevant 
part P\Ut{Kg ® ke)- There is, however, a natural way 
to construct a legitimate dynamical map A t such that 
A-t(KgK s ) is a well defined state. Let a be a system op- 
erator and define the dynamical map At as follows 

(A t («gK^),o) = (K t , (a(g)I E )K t ) 
= (P K t ,(a(E)lE)PoKt) + (PiK t ,(a®l E )PiKt) , (18) 

for all i > and kq = Kg ® Ke- Simple calculation shows 
that 

A t (p) = Z t pZ} + Tr E (Y t [p® n\Y?) , (19) 



where the operator Y t is defined by 

Y t = [ ds e- l(t - s)Hll H 10 (Z s ®I E ) . (20) 
Jo 

Recalling, that A t (p) = Tr E (U t [p®to]U}) with U t = 
e~ %t one has the following relation 

Tv E (U t [p®n}ut) = Z t pZj + Tr E (Y t [p® f2]Y t f ) . (21) 

The pair of equations (fT?]) and ([T9| defines the main re- 
sult of this Letter. Note, that we have a solution for the 
dynamical map A t but the equation of motion (|17[) is not 
for the map itself but for the operator Zt . It should be 
clear that the formalism developed so far is hardly ap- 
plicable in practice due to the complicated nature of the 
unitary operator cxp(— iHut). 

Interaction picture and Born-like approximation. - In 
searching for the reduced dynamics of an open system 
the standard procedure is to pass to the interaction pic- 
ture to remove the free evolution of the system itself [l[ . 
Consider the Hamiltonian of the total composed system 
in the standard form 

H = H + V = H s ®I E +I<g>H E + V . (22) 

Defining V(t) = e tH oty e ~i-H t onc g n( j s ^ ne following for- 
mula for the evolution in the interaction picture 

A t (p) = Z t pZj + Tr E (y t \p®n]Y?) , (23) 

where Y t is defined by 

Y t = [ ds W t , s V w (s)(Z s ®I E ) , (24) 
Jo 

with unitary 

Wt, s = r T exp(-i £ V n (u)du^ , (25) 

and T denotes chronological product. Finally, system 
operators Z t satisfies the following non-local equation 

d t Z t = -iVes{t)Z t - f Mt,sZ s ds , (26) 
Jo 

where V eS (t) = tr E (V(t) • I® fi) and M M = tr B (K M • 
E®f2) with K t , s = V i(t)W t .sVio(s). Note that opera- 
tors Yf, M t and interaction picture counterparts Y t and 
Mt, s contain an infinite number of multi-time correla- 
tions functions which makes the full problem rather un- 
tractable. The advantage of our approach is that it al- 
lows to a series of consistent approximations. By a con- 
sistent we mean an approximation which results in com- 
pletely positive and trace preserving evolution of a den- 
sity matrix, that is, legitimate dynamical map A t . This is 
a notorious problem in the standard Nakajima-Zwanzig 
approach. The simplest approximation consists in ne- 
glecting Vn(i). Roughly speaking Vn(t) is responsible 
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for transitions within irrelevant part of the Hilbcrt space 
Pi("Hs 8> He)- It means that one approximates the evo- 
lution operator Wt, s = Texp(— i J Q Vn (r)dr) by Is ® Is- 
It leads to the second order approximation for the mem- 
ory kernel M M = Tte[I® fiVbi(t)Vio(s) ] and hence it 
is natural to call it Born-like approximation. The big 
advantage of our approach consists in the fact that the 
above approximation leads to the legitimate completely 
positive and trace preserving quantum evolution. 

Example: spin-boson model. - To illustrate our ap- 
proach let us consider well known spin-boson model |l| 
beyond rotating wave approximation (RWA) defined by 

H s = ujqct+o-- , H E = J dku(k)a?(k)a(k) , (27) 

and the interaction term 

V = <r+ ® [a(f) + at (ft)] + a~ ® [a\f) + a(h)] , (28) 

where cv(f) = J dk f(k)a<(k). As usual a ± are standard 
raising and lowering qubit operators. Note that a form- 
factor 'ft' introduces counter-rotating terms. One easily 
computes 

V(t) = e-^°V+ ® [a(f t ) + at(ftt)] + h.c. , (29) 

where time-dependent form- factors read ft(k) = 
e -iu(k)tj!^ an( ^ a smi ii ar formula for h t . Let tip be 
a thermal equilibrium state of the environment at the in- 
verse temperature /?, that is, Qp — e~^ HE /Tr(e~@ HB ). 
In the standard approach one takes zero temperature 
state floo = | vac) (vac |, with |vac) being the vacuum 
state of the field Introducing n-point correlation 

functions (AiA 2 ■ ■ ■ A n )p = Tr (A\A 2 . . . A n Clp) one has 
(a(f))p = ( a Hf))e = and 

(a(f)a(h))p = {a\f)a^h))p = , 

and hence (V(t))p = 0. Let ke be an amplitude of flp, 
1/2 

that is, ke = tlJ U with arbitrary unitary U in He- In 
the Born-like approximation the formula M tjS reduces to 

M t , s =m 1 (t,s)\l}(l\+m 2 (t,s)\2}{2\ , (30) 

with 

mi (t,s) = e-^( t ^[(a(/ t )a t (/ s )) /3 + (a t (ft t )a(/i s ))^], 
m 2 (t,s) = e i ^ t -^[{a{h t )a\h a ))p + {a^f t )a{f a ))p\. 

It proves that within this approximation the dynamics 
of Z t is fully controlled by 2-point correlation functions 
{a(f t )at(f s ))p and (a t (/ t )a(/ s )) /3 and similarly with f t 
replaced by h t . The above 2-point correlation functions 
are time homogeneous due to the invariance of VLp and 
hence rrik(t,s) = mk(t — s). It is, therefore, clear that 
the solution of (|2l)]) has the following form 

Z t = Zl (t)\l){l\+z 2 (t)\2){2\, (31) 



where the complex functions Zk(t) satisfy 

d t Zk{t) = m k {t - s) z k (s)ds , (32) 
Jo 

with 2fc(0) = 1. Interestingly, we have to decoupled 
equations for z k . Having solved for Z t one computes 
Tr s (Y t [p®fi /3 ]Y t t ) , where Y t = f* dsV 10 (s)Z s ® l E . Ob- 
serving that Vio(s) = PiV(s)P = V(s)P due to 
PoV(t)Po = 0, one finds the following Kraus represen- 
tation for the map p Tr E{Y t [p<E>tlp)Y^): 

di(t)a + pa~ + d 2 (t)a~ pa + + a(t)a + pa + + a* {t)a~ pa~ , 

where 

di{t) = [ ds [ dwe- ?; " o(s -" ) L'i(s-u)zi(s)z 1 (w) , 
Jo Jo 

d 2 {t) = f ds f due lulo{s - u) D 2 (s~u)z 2 {s)z* 2 {u) , 
Jo Jo 

a(t) = f ds f dMe~ lt " o(s+u) y4(s-u)zi(s)z*(u) . 
Jo Jo 

with 

D^s-u) = (a(f s )a\f u ))p + (a\h s )a(h u ))p , 
D 2 (s-u) = {a\f s )a(f u ))p + (a(h s )a^h u ))p , 
A(s - u) = {a(fs)a\h u ))p + (J (h s )a(f u )) p . 

Interestingly, the preservation of trace implies that 

d k (t) + \z k (t)\ 2 = 1 , 

for k = 1,2. Actually, this simple condition is hardly 
visible from the definition of dk (t) and the corresponding 
non-local equations for Zk(t). The standard spin-boson 
model [l| is easily recovered if j3 — >• 00, \zi(t)\ = 1 and 
a(t) = 0. In this case the vacuum correlation functions 
simplify to Di(s - u) = (a(ft u )ot(ft s )) vac , D 2 (s - u) = 
(a(/«)a t (/ s ))vac , and A(s - u) = (a(/ u )a t (/i s )) vac due 
to (at {ft)a(f s )) vac = (a^(h t )a(h s )) vac = 0. Now, in the 
white noise limit 

(a(ft s )a t (ft M )) vac = ji5(s-u), 
(a(/ s )a t (/ u ))vac = j 2 S(s-u), (33) 
(a(/ s )a t (ft u )) vac = r)6(s-u), 

with 71,72 > and 77 € C such that 7172 > |?y| 2 , one 
obtains the following formula for the dynamical map 

Pii(t) = e-T lt p 11 + (l-e-^ t )p 22 , 
P22W = e-T 2t p22 + (l-e^ lt )p 11 , 

P12W = e-^^+^ t pi2+?7(t)P2i , (34) 
P21W = e-a^+^l*/^! +T)*{t)p 12 , 

where 7y(t) = , 2 '' . (1 - e -5(T 1+ T 2+4iw °)). It is 
easy to show that in the standard secular approxima- 
tion rj(t) = 0. Note, however, that even in the secular 
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approximation the above dynamics does not correspond 
to a semigroup unless 71 or 72 vanish, that is, the in- 
teraction Hamiltonian contains either resonant or anti- 
resonant terms only. Interestingly, if both 7s are strictly 
positive then asymptotically one has 

pn(i) P22 , P2i{t) -> pu , (35) 

which means that asymptotically the occupations of lev- 
els pn and P22 simply swap. Hence, such evolution does 
not have a proper equilibrium state and displays highly 
non-ergodic behavior (its asymptotic state highly depend 
upon the initial one). 

Conclusions. - We have provided a new approach to 
open systems quantum dynamics. Contrary to the stan- 
dard Nakajima-Zwanzig projection technique in the Ba- 
nach space of density operators it is based on the Fesh- 
bach projection technique in the Hilbert space of ampli- 
tudes. The main advantages of presented approach are 
i) it is based on a much simpler dynamical equation not 
for the dynamical map itself but for the linear operator 
acting in the system Hilbert space. Solving this equation 
one constructs the exact dynamical map. ii) This ap- 
proach, contrary to used so far Hilbert space techniques, 
enables one to work with mixed states of the environ- 
ment, in) Its crucial property is the ability for coher- 
ent approximations. This is a notorious problem with 



Nakajima-Zwanzig equation. 

We have illustrated a simple Born-like approxima- 
tion within the standard spin-boson model beyond RWA. 
Moreover, contrary to the standard approach [l[ we do 
not assume that an environment is in the vacuum state 
and consider the general invariant equilibrium state Qp. 
Let us observe that this approach enables one to perform 
more refine approximations. The general idea consists in 
"reducing" an irrelevant part, that is, the range of Pi. 
One may introduce an approximate projector P* pp such 
that its range defines a proper subspace of the range of 
Pi . For example in the spin-boson model if p* pp projects 
into 1-particlc sector then P? pp V{t)P? pp = and one 
recovers Born-like approximation. Interestingly, Born- 
like approximation docs not result in a Markovian semi- 
group. Actually, even Markovian semigroup on the level 
of Z t leads in general to genuine non-Markovian evolu- 
tion on the level of a dynamical map A t . It shows that 
Markovian evolution corresponds to much more refine ap- 
proximation than simple Born-like one. It is hoped that 
further research reveals new aspects of quantum evolu- 
tion and intricate relations between Markovian and non- 
Markovian dynamics. 

This work was partially supported by the National Sci- 
ence Centre project DEC-2011/03/B/ST2/00136. 
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